Ergodicity breaking in strong and network-forming glassy system 
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The temperature dependence of the non-ergodicity factor of vitreous GeC>2, fq{T), as deduced 
from elastic and quasi-elastic neutron scattering experiments, is analyzed. The data are collected 
in a wide range of temperatures from the glassy phase, up to the glass transition temperature, and 
well above into the undercooled liquid state. Notwithstanding the investigated system is classified 
as prototype of strong glass, it is found that the temperature- and the g-behavior of f q (T) follow 
some of the predictions of Mode Coupling Theory. The experimental data support the hypothesis 
of the existence of an ergodic to non-ergodic transition occurring also in network forming glassy 
systems. 

PACS numbers: 64.70.kj,61.05.F-,61.43.Fs 
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On cooling a liquid below the melting point of its crys- 
tal, its viscosity increases and its molecular motion is 
slowed down. Below the glass transition temperature, 
T g , no more molecular rearrangements typical of a liq- 
uid occur: the system, becoming a glass, is completely 
arrested, at least on the observation time scale. The 
structural (a—) relaxation process describes the dynam- 
ical arrest. Its characteristic time r a exhibits a strong 
temperature dependence, usually mirrored by that of the 
static transport coefficients, such as the shear viscosity, 
r a oc r\. In the glassy state, only secondary relaxation 
processes, f3, remain activated. They are attributed to 
local reorientational motions. 

From the theoretical side, several different micro- 
scopic pictures have been proposed to universally de- 
scribe the dynamical arrest characterizing the glass tran- 
sition [l], 0, H, 0, [|[ . The existence of a trapping effect 
that molecules or group of molecules suffer, due to the 
presence of their neighbors, was formulated both in the 
phenomenological approach of the free volume and in 
the more formal framework of the Mode Coupling The- 
ory (MCT) [| . The MCT is, up to now, the only ap- 
proach which provides a self consistent treatment of the 
particles dynamics. Thanks to the introduction of a non- 
linear microscopic equation of motion, it allows to calcu- 
late the time evolution of & q (t), the normalized density 
autocorrelation function at momentum transfer q. The 
glass-transition problem is approached introducing the 
existence of a dynamic instability at some temperature, 
T c , above T g . The transition from ergodic (high temper- 
ature) to non-ergodic (low temperature) state takes place 
at this temperature [3J. From the microscopic point of 
view, T c corresponds to a critical density from which each 
particle is trapped in the cage made by its neighbors. 
In this picture, only the rattling motion of the particles 
blocked inside the cages is permitted: this motion is the 
microscopic origin of the (3 relaxation Q . 

In recent years, a considerable amount of work has 



been done to extend and solve the MCT equations for 
systems much more complex than simple liquids with 
spherically symmetric interactions for which the theory 
was orig inally formulated 0,0, H ©]. Some universal fea- 
tures, largely independent from the particular system one 
is analyzing, were foreseen. For instance, in a tempera- 
ture region close to T c , a cusp-like behavior of the long 
time limit of the density correlator, & q (t)j—^^f q (T) is 
predicted. In fact, in a liquid f q (T) is equal to zero, 
indicating that at long enough times a given particle 
can leave the transient cage made by its nearest neigh- 
bors. Crossing T c , no more molecular rearrangement oc- 
curs and f q (T) changes discontinuously assuming a finite 
value. In this sense f q {T) is usually considered as an 
indicator of the ergodic to non-ergodic transition Q. 

The Inelastic Neutron Scattering (INS) and Inelas- 
tic X-Ray Scattering (IXS) experiments provide a good 
mean to study the temperature evolution of f q {T) in a 
wide q range. From a experimental point of view, <& q {t) 
and f q (T) are accessible measuring the intermediate scat- 
tering function, F(q,t). In a supercooled liquid, $ g (t) is 
expected to go to zero exhibiting a two-step decay re- 
lated to the (3 and a relaxations respectively. On cooling 
the sample, the structural rearrangements became slower 
and slower, and the decay of $ q (t) related to the a re- 
laxation is shifted to longer and longer times. When the 
temperature is low enough and the a relaxation is frozen, 
the characteristic time of its decay is no more measurable 
and, for long times, Q q (t) assumes a finite value: f q (T). 
Alternatively, the plateau value separating the (3 relax- 
ation from the a one, sited at longer times, can be deter- 
mined in the frequency space by measuring the dynamic 
structure factor, S(q,ui). When the characteristic time 
of the a process is shifted to so long time that l/r a is 
inside the frequency interval defined by the instrumental 
resolution function, the evolution of the a process can 
be directly monitored by measuring the elastic scatter- 
ing intensity. Under this condition, f q (T) is determined 
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FIG. 1: Quasi elastic scattering measurement at q — 3.37A -1 
are reported for selected temperatures, 300 K (circles), 750 
K (stars) and 1000 K (squares). In the inset the normalized 
spectra are reported together with the instrumental resolu- 
tion measured by a vanadium scan (full line). The spectra 
are identical to the instrumental resolution which, due to the 
experimental setup, is somewhat asymmetric. 



by the ratio of the elastic to the total scattered intensity, 
that is f q {T) = S el (q,u = 0) T /S(q) T . 

In the past two decades, the existence of the spe- 
cific theoretical predictions on the T and q behavior of 
fq{T) motivated a large number of experiments and sim- 
ulations aimed at testing the MCT in different systems 
[3 HO, GJ, El El El ■ Recent studies tried also to relate 
the microscopic properties of the systems with the tem- 
perature dependence of their viscosity, when approach- 
ing the glass transition 0, E3, EH O H3| • Up to now, 
the experimental studies performed by INS [lj], E2> EH 
and more recently also by IXS [HI, EI| were carried out 
mainly on glass formers without directional bounds (Van 
der Waals or ionic systems) which show high fragility 
(high increase of viscosity upon cooling [2l(] V. Recently, 
both theoretical and simulation studies [6j, |7j, and ex- 
perimental data [13, H3| , have suggested that the MCT 
can describe the behavior of f q (T) also in systems with 
a spatial organization of the molecules. 

The aim of the present work is to experimentally test 
which MCT predictions, if any, are valid also for strong 
and network forming glassy systems. For this purpose, 
we choose as prototype system the vitreous germania, v- 
GeG^: it is one of the strongest glass mGeo 2 = 20 [25[ and 
it exhibits a rather low glass transition temperature easy 
accessible by the experiments (T g ~ 80QK). The sam- 
ples were obtained by melt-quench process of the Germa- 
nium(IV) oxide crystalline powder purchase by Aldrich 
(purity greater than 99.998%). The system has been in- 
vestigated by neutron scattering experiments. 

Elastic and quasi-elastic measurements were performed 
on the thermal backscattering spectrometer IN13 at the 



Institut Laue-Langevin (ILL, Grenoble, France). The to- 
tal elastic intensity, S e i(q,u> — 0,T), was measured as a 
function of temperature from 20 K to 1100 K. Quasi- 
elastic spectra were collected at significant temperatures 
in the energy range [-0.026 meV to +0.03 meV] with an 
energy resolution Aui = 9/ieV. 

The measurements of the static structure factor, S(q), 
were performed at the diffractometer 7C2 located on the 
hot source of the reactor Orphee of the Laboratoire Leon 
Brillouin (LLB, Saclay, France). The incident wavelength 
A = 0.729A was chosen to access a wide Q range between 
0.4 to 15A" 1 and to minimize the inelasticity corrections. 
The S(q) was investigated in the temperature range from 
room temperature to the liquid phase. 

The two sets of data were reduced using MonteCarlo 
simulation to estimate the multiple scatterin g, t he cell 
contribution and the transmission coefficients [24|. 

In Fig. [TJ where the quasi-elastic spectra are reported 
for selected temperatures, the decreasing intensity under- 
lines the temperature effect of the atoms motion. In the 
inset of Fig. [1] we report the same spectra, normalized at 
the maximum intensity. No changes of the quasi-elastic 
line shape are observed in the whole explored tempera- 
ture range, indicating that there are no detectable ther- 
mically activated processes, at least within the width of 
the experimental resolution. Therefore the a relaxation 
process is always inside the elastic peak and the relation 
AuT a (T) >> 1 is always fulfilled. As mentioned before, 
under this condition the f q (T) can be analyzed accord- 
ing to the previous description. In particular, considering 
that f q (T 0) = 1, f q (T) is obtained as [l3[: 

S e i{q,u = 0) T = f t T )- S ^ T (1) 
S e i(q, lo = 0)t=o 9 S(q)r=o 

in this way each data set is normalized by its low tem- 
perature measurement, thus cancelling all instrumental 
normalization constants. Since our low temperature da- 
tum in the elastic scan is 20 K , we linearly extrapolate 
the S(q,u> — 0)t=20 to S(q,u> — 0)t=o- Room tempera- 
ture is the lowest temperature where the static structure 
factor was measured. Since the diffraction pattern does 
not show any significant change from room temperature 
up to 750K (as shown in the inset of Fig(2]) , we safely ap- 
proximate S(q)T=o to S(q)T=35QK ■ Smooth temperature 
evolution is found in the S(q) for temperatures greater 
than 750 K. In FigJ^l S(q) is presented at two significant 
temperatures (T = 750 K and T = 1000 K) together with 
the ratio, r(g), between the two spectra. r(q) is almost 
constant in the whole range apart at very low q, where 
slight differences are present. 

The obtained f q (T) values are reported in Fig. [3] as 
a function of temperature, for three selected wave vec- 
tor transfers, q. The data follow the same trend, with- 
out any significant deviation, in the whole q range here 
investigated. In the low temperature region, the tem- 
perature dependence of f q {T) is fairly well described by 
a linear behavior. Increasing the temperature a strong 
slope variation is clearly identified. It occurs at the same 
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FIG. 2: Static structure factor S(q) measured at the two rel- 
evant temperatures T = 750 K (line, upper panel) and T = 
1000 K (circles, upper panel). The ratio r(Q) = S ^f=™° 
is also reported (circles, lower panel). The dashed line in the 
lower panel is an horizontal guide to the eye. The errors bars 
are of the same size as the symbols. In the inset the low q part 
of S(q) measured at T=350K (squares) and T=750K (line). 



temperature independently from the analyzed q value. 
In the highest temperature range f g (T) becomes almost 
temperature independent at a value which strongly de- 
pends on q. Interestingly, a very similar temperature de- 
pendence is also found in v-SiC>2 reported in the inset of 
Fig. [3l These measurements were collected using the IN6 
spectrometer installed at ILL, with q ~ 2.6A _1 [2||. As 
previously described, in the MCT, f g (T) is expected to 
display a cusp behavior in the region close to the critical 
temperature, T c , which should be located above T g . The 
theory foresees that close to T c the following relationship 
holds: 



_ r r q + h q {T - T c y 



T <T C 
T>T r 



(2) 



where (3 — 0.5 is expected and f5 and h q are positive 
quantities depending only on g |3|. In order to apply 
this relation the static structure factor S{q) must de- 
pend smoothly on the temperature: this guarantees that 
any change in the temperature behavior is not related 
to a structural rearrangement, but it has a dynamical 
origin. This condition is fulfilled by our data as shown 
in Fig[2] As foreseen by MCT, the presented data sug- 
gest the existence of a break in the trend of f q (T) versus 
temperature . The strong slope variation in the temper- 
ature dependence of f q (T) is located near T g . It is not 
straightforward to identify the temperature of the break 
as T c defined in Eq. 2 since in strong glasses it is ex- 
pected to largely exceed T g [27l ]. However it is evident 
that a transition occurs and the system starts to anoma- 
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FIG. 3: Temperature dependence of f q (T) measured in v- 
Ge02 at three selected wave vector transfers. The dashed 
lines are reported as guides to the eye. 

In the inset f q (T) is reported for the v-Si02 sample measured 
at q = 2.6A -1 using the IN6 spectrometer [26l ]. 



lously change its dynamics. It is quite interesting to ob- 
serve that the same behaviour is observed also in v-Si02 
where T c is expected to be much higher than the temper- 
ature where the break occurs, i.e. 1500 K. Concerning the 
temperature dependence of f q (T), in the cases of these 
two strong glasses the linear trend of the non-ergodicity 
factor is maintained up to temperatures close to T g . We 
have to recall that Eq.(2) is valid when T is close to T c 
(probably a few degrees), while the experimental data 
span a wide temperature range. Indeed, the data are 
perfectly compatible with a power law, with (3 = 0.5, at 
least within some 50 K above T g . This behavior can be 
obtained taking into account that the measured f g (T) is 
the superposition of two contribution: the first related to 
the relaxational process described by Eq.(2), the second 
related to the pure vibrational contribution. The latter, 
calculated b y th e density of vibrational state measured 
at T=300K [23], show an almost linear trend in the low 
temperature range. However, an accurate analysis in the 
region close to T c is probably beyond the present investi- 
gation capability in terms of energy resolution, Q range 
and statistics. 

A further test to ensure that the constant value reached 
by fq(T) in the high temperature range is really fjj of 
eq.2, is represented by its q dependence. According to 
MCT calculations for simple liquids, /° has to follow in 
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FIG. 4: Wave- vector dependence for the fjj parameter (dots), 
compared with the static structure factor, S(q), of v-Ge02 
(circles). 



hase the oscillations of the static structure factor, 5(g) 
This behavior has been experimentally observed in 
different glass-forming systems |10| . The comparison be- 
tween the q behavior of 5(g) and /° obtained by our data 
is reported in Fig. 2J Again, also in the present case the 
ability of /° to oscillate in phase with the main feature 
of 5(g) is confirmed. 

Our results confirm that, also in the analyzed system, 



the MCT is able to qualitatively describe some univer- 
sal aspects of the glass transition. Even if the temper- 
ature, at which the ergodicity breaking occurs, is aston- 
ishing close T g , its existence is observed in v-Ge02 and 
seems to be visible also in v-SiC>2. The existence of the 
ergodicity breaking in network forming glassy systems 
has been previously formulated analyzing the results on 
B2O3 where the f q (T) was investigated by different ex- 
perimental techniques [29| . However, in this latter case, 
fq(T) does not show a well-defined discontinuity in its 
temperature dependence but only a gradual change in a 
temperature range about 200 K wide. Our findings sug- 
gest that MCT predictions could be extended to a wide 
class of materials, such as network forming glassy system. 
This result is somehow unexpected, in fact these systems 
are characterized by a local structure determined by co- 
valent bounds and the significate of caging effect on these 
materials is far from being trivial. 
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